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The spontaneously generated magnetic flux observed by Mannhart et. al. (1996) in asymmetric 
45° grain boundaries in YBa2Cu307-^ was explained by them in terms of a d x 2_ y 2 superconducting 
order parameter and grain boundary faceting. This article argues that twin boundaries which contact 
the grain boundary also play an important role. 
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INTRODUCTION 



Recently Mannhart et. al. jjj and Moler et. al. || have observed spontaneously generated delocalized magnetic 
flux in what they call asymmetric 45° grain boundaries in YBa2Cu307_ a; (YBCO). They have explained the existence 
of this flux in terms of a superconducting order parameter having d x i __ y i symmetry and grain boundaries which are 
faceted. The geometry of their experiment and of their explanation is indicated in Fig. 1 . In the words of Ref . || , "For 
a (i-wave superconductor, a lobe of the order parameter is normal to the boundary on one side, while a node is normal 
on the other. Faceting of the grain boundary produces slight variations of the grain boundary angle, rocking the 
local normal from the positive lobe to the negative lobe, varying the local critical current from positive to negative, 
producing a series of and 7r-junctions and spontaneous magnetization." Ref. pj] notes that the variation of the 
• magnetic flux along the grain boundary is described by the equation 



O ■ X J , 2 = sm[7(y) - a(y)\ (1) 

\o : 

where y is the distance measured along the grain boundary, 7(2/) is the gauge-invariant phase difference across the 
boundary at point y, and Aj is the Josephson penetration depth. The idea that the variations of the grain boundary 
angle cause the local critical current to vary from positive to negative is incorporated into this equation by taking 
a(y) to be zero when the critical current is positive, and 7r where the critical current is negative. The flux in the grain 
boundary, integrated from one end of the grain boundary up to the point y, is given by <&(?/) = < &o7(y)/27r where 
the $0 is the elementary flux quantum. In a realistic model of a grain boundary, Aj would be a function of y (as in 
O , Ref. jjj) but this dependence is neglected in the present article. 

O ■ The purpose of this article is to present an alternative physical model which also leads to spontaneous flux generation 
, in asymmetric 45° twin boundaries, and which could thus also account for the experimentally observed behavior. In 
this alternative model, which also leads to Eq. |j], the changes of a(y) from to tt occur where twin boundaries from 
one side of the grain boundary contact the grain boundary (twin boundaries on the other side of the grain boundary 
do not cause a(y) to change). The details of how this twin-boundary mechanism works depend on the fact that YBCO 
is orthorhombic, and also on the assumption that the twin boundaries in YBCO have odd reflection symmetry (e.g. 
see Refs. ||||). The essential idea is that, given these assumptions, the point y = of Fig. |2|(a) can be viewed as a 
meeting point of three Josephson junctions, an odd number of which must be tt junctions. Under these circumstances, 
it is known from the work of Refs. that the superconducting ground state will have a vortex containing a flux 

■i<I>o at the meeting point of the junctions. 



TWIN BOUNDARY MECHANISM 

Consider the situation shown in Fig. ||(a). On the left side of the grain boundary there are two distinct twins 
separated by a twin boundary and characterized by their Ginzburg-Landau order parameters ipi and ^2 > respectively; 
the superconductor on the right side of the grain boundary is characterized by the order parameter vp . In a Ginzburg- 
Landau model for the grain boundary, the relevant contribution to the superconducting free energy per unit area of 
the grain boundary is 

/> = Cife(^o) (2) 



1 



for points y on the positive y axis, and 



/< = DRe{M* ) (3) 



for points on the negative y axis (e. g. see Rcf. 10|). Now consider a new state obtained by reflecting the sample of 
Fig. ^|(a) in a plane normal to the plane of the figure and containing the twin boundary. Assuming that the YBCO 
superconducting order parameter is of the A\ g (or ux 2 + vy 2 ) orthorhombic symmetry type (e. g. see Refs. |^]J|), the 
new order parameters are related to the old by V>™ eu ' = ip2, V^ 6 ™ = V'l ancl ipo ew = V'o- Since the total free energy is 
unchanged by such a reflection, C — D in Eqs. g and |3[ 

The interpretation of Josephson experiments [pj|-|l"4j on twinned crystals has shown that the superconducting 
state of YBCO has odd symmetry with respect to a reflection in a twin boundary, i. e. that tp2 = —V'l- The argument 
is as follows. In Rcf. |lC| ] it is shown that the superconducting order parameter can either be continous or change 
sign at a twin boundary; for YBCO these two cases correspond to the superconducting state having either even or 
odd symmetry with respect to a reflection in the twin boundary. Now, Josephson experiments on twinned crystals 
determine what is called the macroscopic symmetry of the superconducting state of the twinned crystal (which is 
different from the microscopic symmetry of the superconducting state of a single twin). It can be shown that, if 
the twin boundaries have even reflection symmetry, the macroscopic symmetry of the superconducting state of the 
twinned crystal is that of a tetragonal s-wave superconductor. On the other hand, if the twin boundaries have odd 
reflection symmetry, the macroscopic symmetry is that of a tetragonal d x 2_ y 2 superconductor. Since the Josephson 
experiments on twinned crystals show macroscopic d x 2_ y 2 symmetry, the superconducting twin boundaries 

must have odd reflection symmetry. Further details of this argument can be found in (^||. 

A twin boundary at which the superconducting state has odd reflection symmetry (i.e. for which V'l = —^2) can 
be viewed as a 7r Josephson junction (a 7r junction because of the change of sign of the order parameter there) with 
a very high critical current density. Because of the very high value of the critical current density (which is much 
greater than the current densities occuring in the superconductor) the phase difference of the order parameter across 
the junction is always locked to tt. 

Now note that the point y = in Fig. ||(a) can be viewed as the meeting point of three Josephson junctions. 
Because C = D in Eqs. || and [?] the grain boundary Josephson junctions along y > and y < are either both normal 
junctions or both tt junctions. Also, as discussed above, since V>2 = -V'l along the twin boundary, the twin boundary 
can be viewed as a tt junction. Clearly, an odd number of tt junctions converge at the meeting point. Thus, the stable 
superconducting ground state will contain a vortex with flux i$o at the point y = in agreement with the arguments 
and experimental observations of || and 

As noted in Ref. [fl5|| , the designation of a Josephson junction as a normal junction or a tt junction is a matter 
of convention (the total number of tt junctions which meet at a point is however fixed). In the above, the given 
transformation properties of the order parameter under a reflection in the twin boundary are what imposed the 
particular convention adopted. If the substitution V2 = — V"2 is made, then the twin boundary will become a normal 
junction and the grain boundary junction for y > will become a tt junction. Note that this substitution produces 
an extra minus sign in Eq. |[ it is this additional minus sign that leads to a(y) = tt for y > in Eq. EL while a(y) = 
for y < 0. 

It is also important to note that twin boundaries on the other side of the grain boundary, as in Fig. |^(b), arc not 
expected to lead to variations of a(y). As above, the superconducting free energy per unit area of the grain boundary 
above and below the point zero can be written 

/> = Ai?e(V 2 Vo*) > (4) 

and 



/< = Bi?e(ViVo*) • (5) 

The assumption that the superconductivity of YBCO is predominantly of the type d x 2_ y i implies that A and B have 
opposite signs. In addition, for the reasons stated above, ip2 = -V'l- The combination of these two sign changes means 
that there will be no effect on the net contribution to a(y) when crossing a twin boundary of the type indicated in 
Fig. |(b). 

The arguments just given show that the superconducting grain boundary free energy can be written 




where 



2 



V(y) = cos[7(y) - a(y)] 



(7) 



j c is the critical current density, L is the thickness of the superconducting film in the direction normal to the plane 
of Fig. ||, and c is the speed of light. Also, Aj 2 = <&oc/(87r 2 d7 c ) where d w 2Al is the magnetic thickness of the grain 
boundary, Xl being the London penetration depth. The quantity proportional to A 2 represents the magnetic and 
surface current energy in the grain boundary. || 

The free energy of Eq. (^) has the form of a Lagrangian describing the motion of a particle of mass Aj 2 in a potential 
V(j); here y represents the particle displacement, and y is time. Thus, it is clear that the condition for F to be a 
minimum is the equation of motion (1), where the force on the particle is —dV/dy = sin(7 — a). Since the particle 
velocity, dy/dy, is proportional to the magnetic induction h(y) in the junction, i. e. 

^ = 4^A L %)/$ , (8) 
dy 

it will be continuous across the twin boundaries where a changes by 7r; also, since dy/dy is continuous, 7 will also be 
continuous. This interpretation of Eq. ([!]) as describing the motion of a particle in a potential V(y) allow a simple 
intuitive description of the solutions of Eq. (pi). 



JOSEPHSON VORTICES 



First consider the solution of Eq. |l| for a = 0, i. e. no twin boundaries present. j3j,|4| There will be a solution 
corresponding to a particle starting at the potential maximum (point a of Fig. |^(a)) with zero velocity, accelerating 
down the hill to reach its maximum velocity at point b, and then gradually decelerating as it climbs the hill to come 
to rest again at point c. The corresponding graph of displacement versus time (or 7 versus y) is shown in Fig. ||(b). 
The total change in 7 of 2n gives a flux associated with this vortex of <&o (recall <&(y) = 3>o7(y)/(27r)). Since the 
magnetic induction in the grain boundary is proportional to the particle velocity dy/dy, it has its maximum value at 
the point y = of Fig. |](b). For this example, the solution y(y) and the free energy can be evaluated analytically. 
H The free energy of this vortex (relative to the state with y{y) — constant = 0) is positive, which means that the 
vortex is not stable in the absence of an external magnetic field. 



VORTEX WITH FLUX ±$ TRAPPED BY A TWIN BOUNDARY 



Now consider the case of Fig. ^(a) where the grain boundary is contacted by a single twin boundary at y — 0, and 
there are no other twin boundaries within a distance of many Josephson penetration depths Aj. Here, a — for 
y < and a = tt for y > 0. The relevant potential is shown in Fig. ||(a), and the solution of interest corresponds to 
a particle starting at point a with zero velocity, accelerating down the hill to point b, decelerating up the hill from b 
and coming to rest at c. As noted above, the particle velocity is continuous at point b. The corresponding graph of 
7 as a function of y is shown in Fig. [|(b). The change of 7 by it corresponds to the vortex containing a flux of 
This vortex is confined to be centred on the twin boundary, since moving the vortex away from the twin boundary 
will increase its free energy by an amount which, at large distances, is proportional to the distance it is displaced. It 
is also clear that the state with no vortex (i. e. 7 = constant) is unstable with respect to the formation of the i$o 
vortex state. 

A related problem is that in which there are an odd number of twin boundaries like the one in Fig. |^(a) which 
contact the grain boundary in a region y\ < y < y2, but where there are no changes of a (i. e. no twin boundaries) for 
several Josephson penetration depths for y > ?/2 and y < y\. In this case also, the state with no net flux in the region 
yi < y < t/2 will be unstable with respect to the formation of a state with a flux of ^$0 (states with flux (n + 5)^0) 
n an integer, are also possible). 



PINNING OF A JOSEPHSON VORTEX 



Suppose the grain boundary contains only two twin boundaries of the type shown in Fig. y(a), and that these are 
separated by a distance d <C A,/. The appropriate potential is shown in Fig. || Again the particle begins at a with 
zero velocity, accelerates downhill, passes the points b and c corresponding to the location of the two twin boundaries, 
and eventually comes to rest at d. In passing point 6, the phase a jumps by n; because dy/dy must be continuous at 
b, the kinetic energy is continuous at b; the potential energy of the particle in the region from b to c is therefore taken 
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to be C + cos[7(y) + 7r] where the constant C is chosen so that the potential energy, and hence the total energy, will be 
conserved when the particle passes point b. Because of the time the particle spends in the region from b to c, its free 
energy is lowered relative to what it would be in the absence of twin boundaries by an amount which is approximately 
^>oj c Ld/(nc). This means that the vortex will be pinned by the region between the two twin boundaries. 



RANDOM, CLOSELY-SPACED TWIN BOUNDARIES 



Suppose the separation between twin boundaries is random and typically has a magnitude less than the Josephson 
penetration depth Aj. Mannhart et. al. [0] have discussed a stochastic approach to the estimation of the variations 
of "f(y) in this case, and have shown how a flux of <f>o can be spread out over a large distance (i. e. a distance of 
many times Aj) as is observed in their experiments. This spreading out of the flux can also be understood in terms 
of the dynamical picture discussed above. Consider, for example, a particle beginning with zero velocity at a in the 
potential of Fig. and accelerating downhill to b where it encounters a twin boundary; the particle then climbs from 
6 to c; at c the particle velocity goes to zero (the potential energy at c is the same as at a) and the particle then rolls 
back downhill past b, c and comes to /, where another twin boundary is encountered, and then continues towards g. 
It is seen that the particle can spend a lot of time rolling back and forward in the various potentials it encounters 
each time a twin boundary is crossed. This time is equivalent to distance measured along the twin boundary. This 
shows how a small amount of total flux can be spread out over a large distance. 



CONCLUSION 



The fundamental process which is at the origin of spontaneous flux generation in asymmetric 45° grain boundaries 
in the mechanism advanced in this paper is the stabilization of a vortex with flux at points where a twin 
boundary such as that shown in Fig.0(b) intersects a grain boundary. The existence of such half integral vortices is a 
consequence of the fact that the superconducting state of YBCO has odd symmetry with respect to a reflection in a 
twin boundary, and the observation of such vortices would thus constitute an experimental test of this odd reflection 
symmetry. Furthermore, this conclusion is independent of whether the A\ g orthorhombic-symmetry order parameter 
of YBCO is predominantly d x i_ y i or s-wave. 

In contrast to the above, the point where the twin boundary of Fig. 0(b) intersects the grain boundary should not 
support the existence of a vortex containing a flux of ^$o- in this case, the assumption that the YBCO supercon- 
ductivity is predominantly d a; 2_ y 2-wave has been combined with the assumption that the twin boundaries have odd 
reflection symmetry to arrive at this conclusion. 

The fundamental process which is at the origin of spontaneous flux generation in the mechanism of Mannhart et. 
al. and Moler et. al. 0] is the stabilization of a vortex with a flux of a t the vertex b in Fig. of the intersection 
of two distinct facets of a grain boundary. It is possible that both the faceting mechanism and the twin boundary 
mechanism play a role in grain boundaries studied in Ref. [jj] and [0| . 

Although the materials and sample preparation problems are no doubt formidable, it would be of interest to prepare 
samples having isolated twin boundaries, and an isolated intersection of two facets [such as shown in Figs and §(a) 
and (b)] and to observe directly the presence or absence of vortices with a flux of J<E>o in these idealized cases. The 
superconductor on the right hand side of Figs. and 0(a), and on the left of Fig. 0(b), does not necessarily have to 
be YBCO, but could be an isotropic superconductor such as lead. 

An interesting special case of the faceting mechanism for the stablization of a ^$>o vortex occurs in the corner 
junction such as that used in Ref. |I3[ . To obtain the |$o vortex, the junction would have to be long in comparison 
with the Josephson penetration depth. The actual experimental conditions of Ref. 13 correspond to the short junction 
case. 

Note also that in c-axis tunneling from a twinned crystal of YBCO to Pb (as in Ref. |l(|) the Josephson junctions 
formed by the two different types of twins must be, one a tt junction and the other a normal junction (see Ref. [0]). 
Thus, a necessary condition for spontaneous flux generation is satisfied here also. For example, for a YBCO crystal 
which consists of only two twins, a vortex containing a half a flux quantum would form along the line where the plane 
containing the single twin boundary contacts the junction. The fact that the junctions Jl6[ are "fabricated on a broad 
surface normal to the c axis" however means that the normal and tt junctions are distributed over a two-dimensional 
surface. This is different from the situation of Refs. [0,0] where the normal and tt junctions are distributed along a 
one-dimensional grain boundary. 

The degree of orthorhombicity of the superconductivity of YBCO plays a role in the mechanisms just described. 
For example, for a tetragonal d x 2_ y 2 superconductor, the critical current for an ideal asymmetric 45° twin boundary 
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(such as oo' in Fig. [j]) is zero, and the critical currents on the grain boundary facets ab and be have opposite signs. 
For a superconductor which has some orthorhombicity but is nevertheless approximately d x 2_ y 2, the critical current 
on oo' of Fig. [I] is not zero, and the angle 5 (or 6') will have to exceed some nonzero critical value in order for the 
critical currents of the two facets to have opposite signs. Similarly, some orthorhombicity is necessary for the critical 
current not to be zero on the grain boundaries of Fig. ^. 

One measure of the degree of orthorhombicity of the superconductivity in YBCO, is the ratio of the penetration 
depths, A a /Ab, measured along the crystallographic a and b directions, in untwinned single crystals. p7|^8[ This ratio 
is approximately 1.5 to 2.0 in high quality crystals fl7|| thus indicating a considerable degree of superconducting 
orthorhombicity. The result A a /Ah ~ 1 obtained in an earlier measurement [|l9| has been attributed [ 
presence of disorder in the CuO chains. Since the degree of orthorhombicity appears to depend on crystal quality, 
and since the orientation S of the grain boundary (see Fig. |l|) at which the critical Josephson current vanishes is also 
a measure of the orthorhombicity, the crystal quality is expected to play some role in the detailed understanding of 
spontaneous flux generation in asymmetric 45° grain boundaries. 
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FIG. 1. Two tetragonal d x 2_ y 2 superconductors with their a and b axes in the plane of the figure but oriented at 45° degrees 
relative to each other are shown schematically. The two superconductor are on either side of a grain boundary. The dashed 
line oo' represents an ideal straight asymmetric 45° grain boundary, whereas the solid lines ab and be represent two facets of a 
grain boundary which are rotated by angles of S and 5' relative to the ideal case. 
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FIG. 2. (a) and (b) The two ways in which twin boundaries can end at a grain boundary. The directions of the crystallographic 
a and b axes in the different regions are indicated by the rectangles. 

FIG. 3. (a) The dynamical potential V{"y) appropriate for the description of a Josephson vortex, (b) Qualitative behavior 
of the gauge invariant phase 7 as a function of the distance y along the grain boundary for the Josephson vortex. 

FIG. 4. (a) The dynamical potential V(j) appropriate for the description of a vortex containing a flux |$o- (t>) The gauge 
invariant phase 7 for this vortex. The origin y = is the location of a twin boundary of the type shown in Fig. 0(a). 

FIG. 5. The dynamical potential describing a Josephson vortex pinned by a pair of neighboring twin boundaries. The changes 
in slope of the potential at b and c are due to the twin boundaries. 



FIG. 6. Dynamical potential for which twin boundaries are encountered at points b and /. 
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